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The dynamics near a hyperbolic fixed point in phase space is modelled by an inverted harmonic 
oscillator. We investigate the effect of the classical instability on the open quantum dynamics 
of the oscillator, introduced through the interaction with a thermal bath, using both the survival 
probability function and the rate of von Neumann entropy increase, for large times. In this parameter 
range we prove, using influence functional techniques, that the survival probability function decreases 
exponentially at a rate, n' , depending not only on the measure of instability in the model but also 
on the strength of interaction with the environment. We also show that n' determines the rate of 
the von Neumann entropy increase and that this result is independent of the temperature of the 
environment. This generalises earlier results which are valid in the limit of vanishing dissipation. The 
ON , validity of inferring similar rates of survival probability decrease and entropy increase for quantum 

ON . chaotic systems is also discussed. 

PACS numbers : 03.65.-w, 05.40.+j, 05.45+b 
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I. INTRODUCTION 



Quantum-classical correspondence for the case of classically chaotic systems has been much investigated recently 
(e.g., Refs. Expectation values of corresponding dynamical variables begin to differ [0J^], as do classical and 

On : quantum phase space distributions |^|,[l0|, on extremely short timescales which are typically logarithmic in Planck's 
constant. If these studies are taken at face value, therefore, all chaotic dynamical systems, being fundamentally 
quantum in nature, should either be obeying quantum laws of evolution now or be expected to do so in an extremely 
■ short time. Observations tell us otherwise, however. 

Sarkar and Satchell a decade ago, already pointed out the possible role of environment in the quantum 

0^ . evolution of chaotic systems. Recently ]12j, Zurek and Paz have conjectured an interesting quantitative relation 
between a classical chaotic system and its quantum version which is in contact with a bath. They have considered 
the Wigner representation of the quantum Liouville equation 



h 2n (-i) n d 2n+1 vd 2n+1 w 



o: w = {h, w} PB + ^ 22n{2n + 1}! dx2n+1 dp2n+1 (i) 

J-J ' for a particle in a potential V(x) moving in a 2-dimensional phase space. Clearly the h terms are a singular perturbation 
J> , of the classical Liouville equation, in that the order of the differential equation is changed. For chaotic systems 
derivatives of the Wigner function with respect to momentum become large enough to render the quantum correction 
terms comparable in magnitude to the classical Poisson bracket after the Ehrenfest time, tk cx (1/A) log(l/7i), where 
A is a Lyapunov exponent. Therefore, for t > t% we expect significant differences between the classical and quantum 
descriptions of the same system. 

We will now consider an environment consisting of harmonic oscillators to be coupled to the particle. The Caldcira- 
Leggett model |l3|,|l4| will be used for these oscillators. For the special case of a high-temperature, Ohmic environment 
the right hand side of eqn.(|l|) is modified by the addition of a dissipitive and a decoherence term fL^| . When the 
temperature is high enough and the dissipation sufficiently low the dissipative term may be considered unimportant 
and the decoherence term only survives. The non-unitary Wigner function evolution becomes 

• ^ n 2n {~i) n d 2n+1 vd 2n+1 w d 2 w 

W -{ H > W Ipb + 2^ 22n{2n+iy dx2n+1 dp2n+1 + D ^T- ( 2 ) 

The decoherence term is in the form of a diffusive contribution to the dynamics with diffusion coefficient, D. This 
is vital, since it is the diffusion resulting from the opening of the system which limits the development of the fine 
structure in the momentum direction to a critical momentum scale, o~ c . The timescale on which this process occurs 
is given by jl2| 



i 



where o~ p (0) is the initial width of a Gaussian wavepacket in the momentum direction. Classical behaviour is recovered, 
therefore, provided the environmentally-induced diffusion process can prevent the development of fine-structure, i.e. 
if 

Tc<-T h . (4) 

However, opening a system to a thermal environment has other consequences as well. In particular, the von 
Neumann entropy of the system, given by S(t) = — Trp r (i) In p r (t), where p r (t) is the reduced density matrix of the 
system at time t, will increase; information is lost to the environment and initially pure, superposition states of the 
quantum system become classical mixtures in a very short time. Since our ability to predict accurately the behaviour 
a classical system exposed to an initial perturbation depends very much on the nature of the dynamics, it is natural 
to ask whether this is true for the rate of information loss of its quantum analogue due to a perturbing environment. 

As a first step towards answering this question (see also [[il|[l6| ) Zurek and co-workers have considered the inverted 
harmonic oscillator |l7| of unit mass with Hamiltonian 

(») 

This is intended as a model of instability and, in fact, the dynamical behaviour in phase space is dominated by a 
hyperbolic point at the origin. The unstable and stable directions, and the rate at which initial phase space distibutions 
expand and contract in these directions, respectively, are determined by A. In this sense we call A an "instability 
parameter" analagous to a Lyapunov exponent in a classical chaotic system. Indeed at any point on a trajectory 
the sum of the Lyapunov exponents is zero. For a chaotic trajectory there must be one pair of non-zero Lyapunov 
exponents. 

However, there are a number of reasons why we should question any conclusions drawn as to the implications for 
a real chaotic system based on so simple a model. Firstly, there are no quantum corrections to the Wigner function 
evolution for this quadratic potential. The model does not allow for these influences on the dynamics which, though 
small in the presence of an environment in comparison to the classical terms, nonetheless are generally always present. 
The stable and unstable manifolds associated with all hyperbolic points in Hamiltonian chaotic systems intersect both 



one another and those associated with other hyperbolic points 18 1. In this way homoclinic and heteroclinic points are 
formed. The stable and unstable manifolds of the inverted oscillator intersect only at the hyperbolic origin in phase 
space. Clearly, therefore, the effect that the complicated distribution of homoclinic points might have on the open 
dynamics is not taken into account. Neither, of course, is the effect of heteroclinic points. 

Notwithstanding these objections however, the inverted oscillator remains a tractable model of instability both for a 
closed system and for an open system in the presence of an environment. As such, it deserves attention for the insights 
it might give regarding the qualitative and, maybe, quantitative behaviour of genuine, open quantum analogues of 
classically chaotic systems. 

The entropy production rate has been considered in the limit of high temperature and low dissipation. This entailed 
using the approximate Wigner function evolution in eqn.(|^). Zurek and Paz show (14] that after a time determined by 
both A and the strength of the interaction with the environment the rate of entropy increase approaches a constant, 

S -» A, (6) 

i.e., the quantum entropy production rate is determined, in this approximation, by the classical instability parameter. 
Given that the classical Lyapunov exponent to which A is analogous is equal to the Kolmogorov-Sinai (KS) entropy of 
the system |l9| , this is indeed a remarkable characterisation pofl . It suggests that after a time, a quantum, classically 
chaotic system loses information to the environment at a rate determined entirely by the rate at which the classical 
system loses information as a result of its dynamics, namely the KS entropy pj] . 

In this paper we will examine once more the "toy" model of Zurek and Paz. Apart from entropy production we 
will consider the experimentally relevant survival probability function. We will not be restricted to the assumption 
of low dissipation made by others |]l2| . As we will show, the asymptotic behaviour of both the survival probability 
function and the rate of entropy increase will not be determined by A alone but by A in a specific combination with 
the dissipation parameter (which determines the strength of interaction with the environment). Our approach does 
not use the master equation approach of Zurek and Paz but rather Feynman- Vernon influence functional techniques 
p2| . This allows a straightforward analysis of strong coupling of the system to the environment. 

The remainder of this paper is organised as follows: In Section 2 we define the initial state of both the system and 
the environment, determine the time evolution of the reduced density matrix of the system generally and describe 
also how this allows us to calculate the von Neumann entropy, S(t). We specify the nature of the environment more 
explicitly in Section 3, leaving us in a position to consider finite temperature evolution. In Section 4 we shall define 
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the survival probability function, P(t), calculate it for the inverted oscillator and discuss its significance for quantum 
chaotic systems. In Section 5 we show analytically our generalisation of eqn.(^|) for the finite-temperature case. We 
state our conclusions in Section 6. Finally, the appendix contains the more tedious details of the calculations. 



II. REDUCED DENSITY MATRIX DYNAMICS 



A. Initial state 



We shall consider as our initial state the wavefunction 

ip(xi,0) = (by/ir)~i exp 
for which it is easily verified that 



-{Xj - x Q ) 2 
2b 2 



+ ipoX; 



(7) 



[Xi) = Xq, 

b 2 

(p) = Po, 



(p 2 ) 

(A Xi f(Ap) 2 



^2+p , and 
4 ' 



i.e., a state of minimum uncertainty. However, we shall be concerned with density matrices and their evolution so we 
define the normalized initial density matrix corresponding to ip(xi,0) above by 



= Wexp(— e(x 2 + x'j 2 ) + axi + a*x-), 
where the normalization constant J\f is given by 



AA=(&^)- 1 exp(-^f), 



(8) 



(9) 



with e = (26 2 )~ 1 and a = xo/b 2 + ipo. 

We assume now that the system is put into contact with an environment at time t — 0. We will use as our envi- 
ronmental model a bath of independent harmonic oscillators in thermal equilibrium at inverse temperature 1/ (kgT). 
This allows us to write the initial, uncorrelated state of the combined system in operator form as 



PSe(0) = ps(0)®pE(0), 

where 

{Xi\ps(0)\x'^ = p s {x l ,x' l ,0) 
as defined in eqn.(^), and where pe(0) defines the thermal environment, 

mhujr. 



Pe(0)=Y[ 



1 — exp 



k B T 



\m)(m\ > , 



(10) 
(11) 

(12) 



i.e., in a factorized form because of the choice of non-interacting modes. 
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B. Time evolution propagator 



The initial state being so defined, we now concentrate on the time evolution. We assume the Hamiltonian of the 
combined system to be 



where H$ has been given in eqn. 



Hse = Hs 



2 1 

Ld a 



2 2 

i.e., the Hamiltonian of our chosen bath with canonical commutation relations [q n ,p m ] — iM„ im , and 

Hi = -xc(t) y^g n , 



(13) 



(14) 



(15) 



the Hamiltonian of interaction describing the (possibly time dependent) coupling of the inverted oscillator, through 
its position variable, to the position variable of each of the environmental oscillators. 



The combined system and environment, initially in the pure product state of eqn. (10), will, of course, evolve 
unitarily under the action of the Hamiltonian Hg e of eqn. (O) . We will be interested in the reduced density matrix of 
the system at some later time t > .which we write as p r (t) . To arrive at this expression the path integral method of 
Feynman and Vernon is employed [ |13|j22| . The initially uncorrelated state allows us to calculate the evolution kernel 
for the reduced density matrix in the position basis in the following way: 



p r (xf,x'f,t)= I dxi i dx' i Jr{xf,x / f,t\x i} x' i7 0)ps{xi,x' i ,0) 7 



where 



J r (xf,x'f,t\xi,x' i ,Q) 



Di 



Dx' exp<{ -(S[xj ~ S[x'])}F[x,x'], 



(16) 



(17) 



where S is the action of the system and T is the influence functional. Note that T = 1 in the absence of an 
environment. The propagator, J ri has been calculated explicitly |2^;|2^] by Koks, Matacz and Hu and we quote here 
their result 



5 2 



--(a n A 4 2 + a 12 A f A z + a 22 A 2 f ) 

where we have used the more convenient sum and difference coordinates defined by 

A = x-x', S = (x + x')/2, 



(18) 



(19) 



while 64 and an, a\ 2 and a 22 are time dependent terms defined entirely by the spectral density and temperature 

of the thermal environment. We will define them explicitly below when we specify the nature of the environment 
more precisely. 



C. Reduced density matrix evolution 

The final step in calculating the time evolution of the system is to use eqn.([l^) to calculate the reduced density 
matrix p r (t) in the position basis via eqn. @ and eqn.rtq). One finds, after some lengthy but trivial algebra 

Pr (x f ,x' f ,t) = A= exp (— TiAj - r 2 A/E/ - T 3 Z} 

+ r 5 s / + r 6 A / + r 4 ) , (20) 
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where Af has been defined in eqn.(^) above, 



D = 4ftV + b\ + 8eftcm, 



and where we have made the abbreviations 



«22 1 / /e an\ 2 , hai2b 3 



ft D 1 \2 ft 



61 1 / i> 2 &3& 



2t ^ - ^ 



2ft D \ 2% 



26b 2 . 



2b 2 a 12 e > , 



D 



P *2 , . a 2 ^ 2 / , x , 2anftaa* 
T 4 = a A 7 + -^-(e-yl5) + 



'2aft 2aufta* 



1ft 2a*A 7 



and 



(2ah 2 



2anha* 



V2 2a* A 7 



2a\\ha 
D 



2a\\ha 
D 



but in which we have also defined 



A 5 = --a n n, 



A 7 



(e-A 5 ) 



D 



and 



X2 



!J2 





ib 2 


x x = 


T' 


ibz 


012 


2ft ~ 


ft 


ib 3 


ai2 


2h 


ft 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 

(27) 

(28) 
(29) 
(30) 
(31) 

(32) 



The expression in eqn.(po|) is of a Gaussian form and can be diagonalized pq| . The von Neumann entropy, 

S{t) =-tip r (t) In p r (t), (33) 
can therefore be calculated and written in the form 



S(t) = (polnpo + <?lng) , 

Po 



where po and q are defined by 



and 



using eqns.(E2) and (E4) above. 



2JT, 



Pa 



17+ Via' 



= 1 -Po, 



(34) 



(35) 



(36) 
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III. ENVIRONMENT SPECIFICATION 



A. Generalities 

The influence functional used in the determination of the evolution kernel of eqn. ( |l6| ) is determined entirely by the 
so-called "dissipation" and "noise" kernels of the chosen environment. If we now restrict each oscillator in the bath 
to have equal, unit mass we can write 

/>oo 

fi(s,s') = - I du I(uj, s, s') sin(w(s - s')) (37) 

and 







v(s,s')= I d,LdI(uj : s, s') coth ( ) cos(o;(s — s')) (38) 



for the dissipation and noise kernels, respectively, where 



F[x,x'] =exp|-i^ ds^ ds'A(s) [v(s, s')A(s') + 2i/x(s, s') s ( s ')] f" • i*)) 



Here, /(w, s, s') is called the spectral density of the environment as we have assumed the oscillators to have a continuous 
distribution of frequencies, ui fl3|| . Notice that n(s,s') is independent of the temperature of the environment. 

Restricting the discussion now to the case of the inverted oscillator we can determine the time-dependent quantities 
i>i(t), bsit), bs(t) and 64 (t), as well as an(t), ai2(t) and 022 (t), by 

&l(t) = «2(*), 

6a(*) = «2(0), 
6 3 (t) = ui(i), 

&*(*) = iii (0) (40) 



and 

1 



ds / ds'vi(s)v(s,s')vj(s'). (41) 



1 + % Jo Jo 

The functions it, and which determine these quantities are solutions of the differential equations 



il(s) - X 2 u(s) + 2 ds'fj,(s,s')u{s') = 0, (42) 
Jo 

v(s) - X 2 v(s) -2 1 ds'fi(s, s')v(s') = (43) 



when the boundary conditions 



are imposed. 



m (0) =«i(0) =1, ui(t) = v 1 (t) = 0, 
ua(0) =«a(0) = 0, ua(t) = = 1 



B. Calculation of foi(i) 

It is clear from eqn. (^p|) and eqn.(^) above that the quantities b\(t) to 64 (t) depend only the dissipation kernel, 
/x(s, s'), and not on the temperature of the environment. We will concentrate on spectral densities of the form 

a, s') = ^ exp (-—) c (s)c(s'), (44) 



7T V W. 



G 



i.e., densities of an Ohmic type with an upper or cutoff frequency ui c [ fL3| . Using eqn.ji^) in eqn.(^) we find 

M(*V) -> 2 lQ c{s)c{s')5'{s - s') (45) 

as lo c — > oo. A further restriction to the case of constant coupling constants, i.e., c(t) = 1 Vi, enables us to write for 
the inverted oscillator 

b\(t) = (— 70 + Kcoth nt), 
bi(t) = (—70 — recoth/rf), 

Kexp( 7o i) 

&2 (* J = . r 7 and 

smh ra 

The shorthand definition 



«=V A2+ To (47) 

has been used here and we will see that it is an important quantity in the sections below. The asymptotic behaviour 
of bi, i = 1, . . . , 4, in eqn.(|4^) for large nt can easily be determined and each is given by one of 

&!(*)-► (-70 +«), ( 48 ) 
b 2 (t) -> 2 K exp(-( K - 7o )t), (49) 

63 (t) -> 2Kexp(-(K + 7o)t), or (50) 
6 4 (t) -» (-70 - «)• (51) 

Finally, t>i(i) and 1^2 (i) are required for the calculation of each eijj in eqn.(|4l|). Again for the inverted harmonic 
oscillator they are given by 

sinh(/c(t- s))exp(7 s) 

= ■ w 7a (52) 

sinh(Ki) 

and 

sinh(Ac(s))cxp(7o(s-t)) 
W2(S) = sim^) ' (53) 

for s € [0,4], where t is the time at which we wish to calculate the reduced density matrix. 

C. Calculation of the finite temperature (14(f) coefficients 

To examine the asymptotic rate of entropy production of the inverted oscillator, in an environment at a finite 
temperature, we must, of course, calculate the approprate ay coefficients. This will allow to analytically derive the 
asymptotic rate of von Neumann entropy increase. The lengthy details are relegated to Appendix A. Writing A := X/k 
we find eventually 

70 r, s ( dx{n)h(z) 



fll2(z ' 0) = ,m ■ u ^2 E /(") ( d i(»)/s(*) - d 2 (n)U(z)) , (55) 
2(Asmhz) n=(j 

-ioz ™ 

° 22(Z ' 0) = ^ ■ u ^2 E ( d * (")/«(*) + d 2(n)/ 6 (z)) , (56) 
(2Asmhz) n=Q 

where the functions f%(z), i = 1, . . . , 6 are defined in Appendix A, as are d\{n) and ds^ri). Also given in Appendix A 
is the asymptotic behaviour of each (finite-temperature) Oy . 
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IV. SURVIVAL PROBABILITY FUNCTION 



A. Definition and context 



In this section we will consider the survival probability function, P(t), defined by 

P(t) = Tr[p(0)p(i)], (57) 

for a system initially described by the density matrix p(0). In particular, we will be interested in the asymptotic 
behaviour. 

The function P(t) has been examined before in the context of quantum chaos. Tameshtit and Sipe p7j ] have 
considered its behaviour in time for both regular and chaotic systems coupled in a non-demolition fashion to a high- 
temperature thermal reservoir. Using a master equation approach and random matrix theory they show that when 
suitably averaged, the behaviour of P(t) displays substantial differences depending on the nature of the underlying 
dynamics. This function has also been examined theoretically for pure states |28| , p9|] , even for the inverted oscillator 
model p7| . However, the fact that it is accessible experimentally pp| ] makes it a particularly important quantity to 
examine. 

For pure initial states, p(0) — \<p(0))((j)(0)\, and unitary evolution, \<fi(t)) — U (t)\4>(0)) , we have 

p(t) = \(m\m\ 2 , (58) 

i.e., the probability of the system being in the initial state at a later time t. P(t), as defined in eqn. (|57|) above, is a 
generalisation of this applicable to systems which may not be initially pure and/or for which the evolution in time is 
non-unitary. 



B. Survival probability for the open inverted oscillator 



We will now examine the effect of the thermal bath on the survival probability function of the inverted harmonic 
oscillator. For simplicity we will consider our initial state to be centered at the origin of phase space, i.e., we take 
our initial wavefunction to be given by eqn.(Q) with xq — po = 0. The initial density matrix then given by eqn.(|^) 
with a = and Af — (by/n)^ 1 . These initial conditions simplify the form of the reduced density matrix considerably, 
leading to eqn.(p0|) with r4 = T$ = Tq = 0. The survival probablity function can now be simply calculated according 
to eqn.(p7|) above. We find 



P(t) = Tr[p s (0)p r (t)} 



dx 



dx'(x\p s (0)\x')(x'\p r (t)\x) 

1 

r 3 



6tt 2 {(e + 2T 1 )(e + T 3 /2)-T 2 2 /4 



(59) 



where we have yet to specify the magnitude of the parameters of the environment. For the finite temperature case, 
we can easily determine the asymptotic behaviour. Using the asymptotic behaviour of each bi(t) (given by eqns.(|48|) 
- (j5l|)), the asymptotic behaviour of each Oij(t, 0) (given in eqns.(|90|), (91) and (92)), and the definitions of each Ti 
(eqns.(p2|), (pSj) and (p4|)) we easily see that, for a fixed temperature 



{(e + 2r 1 )(e + r 3 /2)-rl/4} 3 ->Ci, 



(60) 



for large times, where C\ is a (temperature-dependent) constant. Consequently, the large time behaviour of P(t) is 
determined entirely by that of T 3 . As t(z) — * 00 we see from eqns.(51), ( |90"| ) and ( pl| ) that D goes to the (temperature- 
dependent) constant 



D a 



sijin 



= 4h 2 e 2 + (70 + k) 2 - (4e?i7o) coth 



(70 - n)n 

2k B T 



(61) 



Now, using its definition in eqn.(|2j 
finally, 



), the asymptotic behaviour of 62 in eqn. (Eg) and that of D in cqn.(|6l|) we find, 
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P(i)*=°°C7 2 exp(-( K -7o)i), (62) 

where C'2 is another temperature-dependent constant. We have checked the accurac y of this result numerically by 
calculating \og[P(t — 1)/P(t)\, which, if eqn.(|62|) is to be believed, has K — 70 = k' = v/A 2 + 7q — To a s its asymptotic 
value. Excellent agreement was found. This is a generalisation of the result of Heller [|6| to the open quantum inverted 
oscillator and reduces to it as 70 — ► 0, i.e. as the system interacts more weakly with the environment, as required. It 
would be natural to conjecture in the spirit of Zurek and Paz fl2]| that the behaviour in eqn.(|6^) might be expected 
to hold for quantisations of classically chaotic systems in interaction with an environment. 

V. ASYMPTOTIC RATE OF ENTROPY INCREASE - ANALYTICAL RESULTS 

We are now in a position to derive the rate at which the von Neumann entropy will increase when z — nt is large, 
i.e., at long times and/or when the Lyapunov exponent is large, for the finite temperature case. 
We can rewrite eqn. (|35|) as 

2a . . 

Po = —, (63) 



where we have made the abbreviation a := yT^/Ti. Also, eqn.(|34|) gives 

dS d ( I 

-tt = -77 i (Po Inpo + (1 - Po) m(l - Po)) 

dt dt I po 

= ••• =^ln(l- Po ), (64) 
Po 

with the overdot denoting a derivative with respect to t. Combining eqns.(|63|) and ( |64| ) we find, for all t, 

? = Aln( 1 ^V (65) 
dt 2a 2 \ 1 + a J y J 

Therefore, both the entropy, S(t), and its rate of change, dS/dt, are determined entirely by the time-dependent 
coefficients Ti and which arise in the expression for the reduced density matrix, eqn. (|2C|) . (Note that these are 
both independent of the centre of the initial minimum uncertainty wavepacket, (xq,pq).) 

The asymptotic expression for D given by eqn.(|6l|), along with that of eqn. (ff9|) , gives the asymptotic behaviour of 
T 3 defined in eqn.(gj): 

8en 2 

lim r 3 = — exp (-2(k - 70)*) . (66) 

nt^oc J-^asym 



Inspection of the asymtotics of the various environmental terms used to define Ti in eqn.(|22|), in particular eqns.(|92| 
©,©,©, © and ©, give 



lim ri -«2 2 (t->oo,0) 



= ^coth ( 

2h V 2k B T 

i.e., also a (temperature-dependent) constant. 

The asymptotic behaviour of a can now be determined using eqns.(p6p and (k 



l„o«h(<2^j, (67) 



lim a = £ exp (— (k — Jo)t) , (68) 

Kt — >00 

where ^ is a temperature-dependent, but time- independent constant. Clearly, a — > as t — > 00, which means that 
Po — > too. Considering eqn. (|6~i|) as t — > 00 we find 







— d 



a(l + a)' 

where we have used eqn.(B3h in the last step. Finally, if we use the asymptotic expression for a, eqn.(pq), we find 



(69) 



dS nt—^oo 

Tt = K ~ 7 ° 



(70) 



This result gives the asymptotic rate of entropy increase in situations where energy dissipation is important and 
cannot be neglected. As 70 determines the strength of the coupling to the environment we can now see that the rate 
at which an unstable, possibly chaotic system will lose information to the environment does depend on the coupling 
strength [p6[ . Note, however, that this asymptotic rate does reduce to the asymptotic rate found previously in the 
weakly coupled regime 



\A 2 + 7o - 7o 

as required. 



dS 
~dt 

« A, when 70 < A, (71) 



VI. CONCLUSION 

In this paper we have used the inverted harmonic oscillator to model instability in open quantum systems. We have 
found that both the survival probability function and the von Neumann entropy increase depend, for large times, on 
the degree of instability in the system and on the strength of interaction with the environment in a simple way. The 
purpose of studying such an elementary system is to build up some degree of intuition as to the behaviour of quantum 
chaotic systems coupled to an environment. A priori the claims of applicability of an inverted oscillator to modelling 
a chaotic system should be treated with caution. We will now discuss to what degree the results for the oscillator can 
serve as a guide to actual quantum behaviour in chaotic systems. 

With regard to the survival probability function prediction we mention the study by D'Ariano et al. jn) of classical 
and quantum structures in the kicked-top model. No environment was included in their model but they have shown 
that the survival probability function (autocorrelation function in their paper) provides an excellent means with which 
to compare classical and quantum invariant structures (periodic points), at least for predominantly regular kicked 
tops. This situation changes when most of the tori have been destroyed. We are currently examining the influence a 
perturbed evolution has on this correspondence. Practically, the survival probability function can be obtained from 
an experimental spectrum (see Heller in |6), also f30| ) and is a useful concept in studies of molecular spectra. The 
longest events in the time domain determine the broadening of the peaks in the spectrum for an unstable periodic 



orbit, with the width being given by the classical Lyapunov exponent. Our result, eqn.(62), suggests that this may 
change in the presence of an environment due to 70. 

For the entropy production prediction we mention two studies of the quantum kicked rotor model jl]|2j|7| evolving 
as it interacts with a perturbing environment |20|]32[ | . In the first po| , the constant von Neumann entropy production 
rate was calculated for various values of the nonlinearity parameter, K, of the system. The Lyapunov exponent, A 
can can be approximated by A w log(/^T/2), for large K and a linear dependence of the constant entropy production 
rate on the Lyapunov exponent was found, i.e. 

S = aX + b, (72) 

where a and b are constants determined by the environment. Moreover, this behaviour is seen in the mixed phase 
spaces resulting from low values of K provided one considers only local Lyapunov exponents to be relevant. In the 
second study p2| an initial state is perturbed unitarily as it evolves in time with the possible perturbations at each 
time step being taken from an ensemble. Averaging over the possible realisations requires that the state of the system 
be described by a density matrix as it is mixed. In this way the entropy can increase. It is shown that the eventually 
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constant rate at which an initial coherent state produces entropy depends to a great degree upon where in phase space 
the centre of the wavepacket is situated: the more chaotic the area of initial localisation - as quantified by a locally 
averaged Lyapunov exponent - the larger the constant entropy production rate. 

We conclude, therefore, that there is indeed some value in using the toy model considered in this paper and others 
Jl2| as a model of instability in open quantum systems. The conjectures which follow from it should be tested in more 
systems which are classically chaotic. 
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APPENDIX A 



In this appendix we will give the explicit definitions of the functions /j, 1, . . . , 6, di(n) and d-zin) used to write the 
finite-temperature expressions for each <Zjj in eqns.(|54|) - (|5^). We will also derive their large-time limits, eqns.(|90|) - 

Once again we will choose constant coupling constants in the ohmic spectral density of eqn.(f44]) and we will also 
assume a large, but finite cj c . This allows us to write 

7(w,s, s)w . (73) 

To calculate the finite temperature noise kernel of eqn. d38| ) we will now formally expand the coth function in a power 
series [ f?3| and integrate term by term. Formally then 

E 9r (*£?) ~ f d " " 2 " cos " (s ~"' ) ' (74) 

Bin are Bernoulli numbers. Defining 

f(r,\ ■= - 

(2n)! \2k B T 
we can rewrite this formal expansion as 

oo 

v{s,s') = 2 1() Y J f^ (2n) (s-s l ), (76) 

i.e. as an infinite sum of derivatives of the (5-function. This expression can now be used in eqn. @, with eqns.(p 
and (|53]), to calculate each a.y. Writing A := X/k wc find 



where 



and 



an(z, 0) = ° ^ £ /(„) + da(n)Mz) , (77) 

(2Asmhz) 2 ^ \ 7o / 

ai2(z ' 0) = ,a ■ u ^2 E (iW^W - d2(n)h{z)) , (78) 
2(Asmhz) n=0 

fl22 ( Z ' °) = Mi ■ : ^ 9 E /(») WW/^) + d 2{n)h{z)) , (79) 
(2Asmhz) 2 n=0 



{ , } (n) := — ((7o - l) 2n ± (7o + I) 2 ") , (80) 
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fi(z) := To 2 cosh2z + 7osinh2z - exp(27 z) + 1 - To 2 , (81) 

h{z) := To sinh2z + cosh2z — exp(27bz), (82) 

f 3 (z) := coshz(exp(2foz) - 1) - To sinh z(exp(27 2:) + 1), (83) 
fi{z) := To coshz(l — exp(2"fo;z)) + sinh z(exp(2To^) — 1) 

+2To 2 sinhz, (84) 

fs(z) := 1 — A 2 exp(27o^) — To ex P(2To- 2 ) (To cos h2z — sinh2z) , (85) 

/ 6 (z) := To (exp(2"f ;z)(To sinh2z — cosh2z) + 1) . (86) 



Asymptotics of a; 

It is easy to show that 



fl{ z ) z^oo _ / „ 2 , » \ / fT\ 

^ 2 (70 + To) (87) 



sinh 2 z 



and 

—2 1. 2 (to + 1) . 

smn z 

These expressions allow us to determine the asymptotic behaviour of au(z, 0) defined in eqn.(|54|). We find 
Iim a u M) = E /(») f ^ + 7o)] + [2(7o + 1)] 



To (To + 1) 



2A 2 n=0 



£/(n)(fo-l) 2 " K 2 ". (89) 



One next uses the definition of f(n) in eqn. (|75j) and the formal power series expansion of the coth function derive the 
desired expression given in eqn.(|9(]) below. To derive eqns.(pl]) and ( p2| ) we proceed along similar lines: 

ail(z , 0) ^_%: coth f(2^_^V (90) 



2 V 2k B T 

a\2 



( Zl Q) _«-dHb). ^^f^^^ (91) 

a22 (z,0)->^cothf%±#V (92) 



2 V 2fc s T 
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